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We investigate the self-locomotion of an elongated microswim- 
mer by virtue of the unidirectional tangential surface treadmilling. 
We show that the propulsion could be almost frictionless, as the 
microswimmer is propelled forward with the speed of the back- 
ward surface motion, i.e. it moves throughout an almost quies- 
cent fluid. We investigate this swimming technique using the spe- 
cial spheroidal coordinates and also find an explicit closed-form 
optimal solution for a two-dimensional treadmiler via complex- 
variable techniques. 

self-locomotion | creeping flow | motility | propulsion efficiency 

Tiny swimmers, be they micro-organisms or microbots, live in a 
world dominated by friction 1 1|. In this world, technically, the 
world of low Reynolds numbers, motion is associated with energy 
dissipation. In the absence of external energy supply objects rapidly 
come to rest |2|. It is both conceptually interesting, and technologi- 
cally important, to try and understand what classes of strategies lead 
to effective swimming in a setting dominated by dissipation. A par- 
ticularly promising class of strategies is where the motion is, in a 
sense, only apparent; where a shape moves with little or no motion of 
material particles. 

The wheel is the mechanical application of this strategy and it is 
instructive to examine it from this point of view. The (unholonomic) 
constraint of rolling without slipping comes about because of the large 
friction between the wheel and the surface supporting it. Neverthe- 
less, a wheel can roll with little or no dissipation of energy. One way 
to view this is to note that the motion of the point of contact is only 
apparent. The point of contact moves, even though the part of the 
wheel in contact with the surface, is momentarily at rest. 

An example closer to the world of low Reynolds numbers, is the 
actin-based propulsion of the leading edge of motile cells [31, intra- 
cellular bacterial pathogens 1 4 1 and biomimetic cargoes 1 5 , 6 . 7 1 . The 
actin filaments assemble themselves from the ambient solution at the 
front end and disassemble themselves at the rear end. Here again, it 
is only shape that is moving and in principle at least the energy in- 
vested at the front end can be recovered at the rear end. (There are 
thermodynamics and entropic issues that we shall not consider here.) 
Apparent motions seem to be good at fighting dissipation (8). 

Here we shall focus on a closely related mode of locomotion: sur- 
face treadmilling. In surface treadmilling the swimmer moves without 
a change of shape, by a tangential surface motion. Surface is gener- 
ated at the front end and is consumed at the rear encfl In contrast to 
actin and microtubules, the surface treadmilling does not rely on the 
exchange of material with the ambient fluid. (The swimmer needs, 
of course, an inner mechanism to transfer material from its rear to its 
front). It is intuitively clear that a needle shaped swimmer undergoing 
treadmilling can move with very little dissipation because the ambi- 
ent fluid is almost quiescent and there is almost no relative motion 
between the surface of the swimmer and the fluid. One can not make 
treadmilling completely non-dissipative because there is always some 
remanent dissipation associated with the motion of the front and rear 
ends. The main question that we shall address here is how can one 
quantitatively estimate this remanent dissipation. 

Let us first consider simple qualitative estimates of the power 
dissipated by viscosity in treadmilling. Consider, a rod-like slender 



swimmer of length £, thickness d with rounded caps undergoing tread- 
milling at velocity U . It is reasonable to assume (and the analysis in 
the following section can be used to justify) that all the dissipation is 
associated with the rounded ends. Hence, by dimensional analysis, 
the power dissipated in treadmilling is of the order of [idU 2 where fi 
is the viscosity coefficient. Let us compare this with the power needed 
to drag the "frozen" treadmiler. By Cox slender body theory 1 12 1 the 
power needed to drag the tube with velocity U is ~ fi£U 2 / log (2£/d). 
Hence the ratio of power invested in dragging and swimming scales 
like (elog l/e) _1 and can be made arbitrarily large. Here e = d/£ 
is the aspect ratio of the swimmer. 

One can now ask if there are slender treadmilers that are arbi- 
trarily better than the slender rod-like treadmiler above? Consider 
now an elongated ellipsoidal microswimmer whose surface is given 
by z 2 /b 2 +r 2 /a 2 = 1 where r 2 = x 2 + y 2 . Let us assume again, that 
the viscous dissipation is a result of a tip propulsion, and estimate the 
position of the tip from the condition |dr/dz| = 1. It can be readily 
demonstrated that in the case of e = a/b <C 1 the tip is located at a 
distance of 6(1 — |e 2 ) from the center and its typical width scales as 
ae. Therefore, applying the same arguments as before, the dissipation 
rate should then scale as V ~ fj,aeU 2 and the ratio of power expanded 
in dragging and swimming becomes b>iU / lo g ( 2ft / a ) ^ 1 

We shall see, by a more accurate analysis, that for prolate spheroid 
the ratio of power in treadmilling to dragging is of the order of 
(e log e)~ 2 . In the following sections we shall analyze two models of 
treadmillers in 3 and 2 dimensions, respectively. 

The theoretical framework. We model the micro-swimmer as a 
prolate spheroid swimming in an unbounded fluid by continuous tan- 
gential surface motion. The Cartesian-coordinate system (x\,X2, £3) 
is fixed with the center O of the spheroid. A modified orthogo- 
nal prolate spheroidal coordinate system (r, ip) with unit vectors 
(e T ,e^,e v ) is defined via the relations x\ — c{t 2 — 1} 1//2 {1 — 
C 2 } 1/2 cos ip, x 2 = c{r 2 - 1} 1/2 {1 - < 2 } 1/2 sin ip and x 3 = cr(, 
where — 1 < ( < l<r< 00, < ip < 2n and c is the semi- 
focal distance |9|. The coordinate surfaces r = r a = const are a 
family of confocal spheroids, x\jb 2 + {x\ + x^/a 2 = 1, centered 
at the origin with major and minor semi-axis given by b = cr a and 
a = c{t 2 — - l} 1 / 2 , respectively. We assume that a steady axisymmet- 
ric flow has been established around the micro-swimmer as a result of 
the tangential surface treadmilling with a uniform far-field velocity U 
(equal to the laboratory frame propulsion speed) in the negative in- 
direction. The low-Re incompressible flow is governed by the Stokes 
and continuity equations, 

Av = n grad p , div v — , [ 1 ] 

respectively, accompanied by the boundary condition at the swimmer 
surface r = r a 

v = u(C) e c . [2] 
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1 Alternatively, we may think of a a slender microbot that is topological^ equivalent to a toroidal 
swimmer proposed by Purcell Q, i.e. the surface is not created or destroyed, but rather under- 
goes a continuous tank-treading motion 
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Since the flow is axisymmetric we introduce the scalar stream-function 
$ (unique up to an additive constant) that satisfies the continuity equa- 
tion 

v — v T e T + v^e^ — curl ( -i— >I> e v ] . [3] 

The velocity components are readily obtained from as 

1 9* 



1 9* 

v T = ^ - — ,«c 



■ — — — 7TT 

H C H V dC 



H Q H V dr ' 



where the symbols H stand for the appropriate Lame metric coeffi- 
cientstf T = c(T 2 -C 2 )i(T 2 -l)-i,H ( = c(r 2 -C 2 )i(l-C 2 )-*, 
= c(r 2 -l)3 (l-( 2 )2. The vorticity field can be obtained from 
E]]as 

u = curlw = -^-E^Ve^ , [4] 



J-£ 2 *e ¥ 



where the operator i? is given by 



E 



c 2( r 2 _ £2) 



r 2 n 92 



Following the standard procedure, the pressure is eliminated from 
the Stokes equation by applying the curl operator to both sides, with 
conjunction with (01] this yields the equation E 4 *? = for the stream- 
function. The boundary conditions []2]] at the microswimmer surface 
t = r a in terms of the stream-function become 



* = 0. 



9 T * 



-<? K-C 2 } ! {i-C 2 } 5 «(C) 



[5] 



and the conditions at infinity (r — > oo) are v T 

The solution for ^ that is regular on the axis and at infinity, and 
also even in £ can be derived from a general semiseparable solution 

dED, 

oo 

* = -2c 2 <7 G 2 (t)G 2 (0 + {A m H m (r)G m (0+ 

m=2,4,... 

Cmtt m (TX)}, [6] 

where Q, m (r, C) is a solution of E 4 ^ = composed from spheroidal 
harmonics that decay at infinity, and G m and H m are the Gegenbauer 
functions of the first and the second kind, respectively. The coeffi- 
cients A m in can be expressed in terms of C m and U via the use 
of the boundary condition = at r = r a . Substituting ^\ into 
[|5ll we arrive after some algebra at the tridiagonal infinite system of 
equations for U and the coefficients C m , 



c(-)/~t I p(°)/~! 



p(+) n 

C-m l^m+2 



ra > 2 [7] 



Here Co 



-c 2 U, E, 



(o,±) 



are known functions of r a , and 



b m — — m(m— l)(2m— 1) 



C 2 



U (C)G m (C)dC [8] 



Regularity of *!/ implies that the admissible solution of []7]] should sat- 
isfy GmHmija) — ► as m — ► oo while the exponentially growing 
solution with C m ~ 0((r a + V T 'i ~ l) 2m ) should be discarded. 

The viscous drag force exerted on the prolate spheroid (in the X3- 
direction) is solely determined by the C2-term in fl6]] corresponding 
to a monopole (Stokeslet) velocity term decaying like 1/r far from 
the particle, F — — (47r/u/c) C2 1111 . Either F or U can be specified 
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Fig. 1. The propulsion velocity of the 'cigar-shaped' microswimmer vs. the 
elongation. 



in addition to the surface velocity, u(C). In the swimming problem 
F = C2 = 0, and u(C) determines the propulsion velocity U . 

The problem of the "frozen" spheroid {i.e. u(Q = 0) in the uni- 
form ambient flow — Ue.z corresponds to substituting b m = 0, ra > 2 
in EJ. In this case the equations for C m can readily be solved yielding 
the well-known result for 9 and the drag force [ TTI . 



F = 



8-kc/j.U 



(l+T,2)cOth 1 T a -T a 



[9] 



Propulsion velocity. In order to determine the velocity of propul- 
sion of the microswimmer freely suspended in the viscous fluid, one 
must solve Eqs. []7]] with C2 = as the particle is force (and torque) 
free. Let us consider the following velocity distribution at the bound- 
ary 

u(C) = -2r a u s (r 2 - C, 2 )~ h (l - CT* Ga(0 > [10] 

where u s is a typical velocity of surface treadmilling and G2 (C) = 
|(1 — C, 2 ). One may verify that for a sphere (c/a — > 0) u(£) = 
w s sinf9, while for an elongated swimmer it(f) ~ u s almost every- 
where except the near vicinity of the poles £ = ±1. More generally 
it can be readily shown that the solution satisfying flU, 111011 is given 
by flU] with Vra , C m = 0; m > 4, A m = and 

A2 = — 2c 2 u s r a (—1 + t 2 ), 

U = U 3 {t 2 - T a (-1 + T a 2 ) COth" 1 T a ) . [11] 

Note that U can be related to the surface motior0 via the use of the 
Lorentz reciprocal theorem [12], 



F U = - / (a ■ n) ■ udS , 



[12] 



where (2, <?) is the velocity and stress field corresponding to transla- 
tion of the same shaped object when acted upon by an external force 
F. For purely tangential surface motion considered in this work we 
have (a ■ n) ■ u = ovc u {0> where 



Vt 

~H~ T 



H T dr 



We calculate the local tangential stress component a T ( from the so- 
lution corresponding to streaming past a rigid prolate spheroid, while 



-In the laboratory frame the velocity at the surface is a superposition of the translattonal velocity 
U and purely tangential motions u 
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r is given by & Substitution into [Jill] with AS = H V H ( dip d( 
yields after some algebra 



U = ~T 



c 2 



u(C) dC, 



[13] 



which holds for an arbitrary tangential boundary velocity u(£). In 
the special case of ) given by II 1 Oil it can be readily demonstrated 
that integration yields the propulsion speed 111 111 . Also, 111311 actu- 
ally solves the infinite tridiagonal system |17|| , since knowing U, F 
(i.e. Co, C2) one can iteratively obtain all the other C m 's by direct 
substitution. The scaled swimming speed of the microswimmer is 
depicted in Figure [T] as a function of the scaled elongation. The val- 
ues of the propulsion velocity corresponding to a spherical swimmer 
(c = 0) and a slender swimmer (c ^> a) can be determined via II12II 
without invoking special spheroidal coordinates. For a sphere, the 
local traction force <r • n = — |f 17 and F = —Q-KfiaU and thus the 
self-propulsion velocity can be found as 1131 



U = 



Aira 2 



u ■ e dS, 



where e is the unit vector in the direction of locomotion. Sub- 
stituting u = u s sin 9 eg and e = e r cos 9 — eg sin 8, where 

9 is the spherical angle measured with respect to e, we arrive at 

2 
3 



U — § Jq Us sin 3 9 dO — | u s in agreement with the result shown 



in FigureQ] 

The drag force exerted on the rod-like microswimmer upon trans- 
lation along its major axis with velocity U ' is given by F w 
— 47T/iaL?'"/(£ log 1/e) and the local friction force is <? ■ n w 

-^"/(alogl/e), where e = a/6 = [1 + (c/a) 2 ]- 1/2 < 1 is 
the aspect ratio. Thus, from 111211 follows 



V 



4:Tvab 



u ■ e dS. 



For the 'needle-shaped' microswimmer the surface velocity u — 
u(Qe( ~ —Us e over almost the whole surface, it follows that 
U ~ u s . As seen in Figure Q] the propulsion velocity U/u 3 — > 1 
as c/a grows and equals to 0.95 already at c ~ 5.3a. As intuitively 
expected, the micro-swimmer is self-propelled forward with the ve- 
locity of the surface treadmilling , while the boundary velocity in the 
laboratory frame is (almost) zero. 

Swimming efficiency. Since the fluid around the elongated mi- 
croswimmer propelled by continuous surface treadmilling is almost 
quiescent, except for the near vicinity of the poles, it is natural to 
expect low viscous dissipation and high hydrodynamic swimming ef- 
ficiency. Several definitions of hydrodynamic efficiency have been 
proposed G2][T4][T9) here we follow the definition 5 = F- U/V, 
where V is the energy dissipated in swimming with velocity U, and 
the expression in the numerator is the work expanded by dragging the 
"frozen" swimmer at velocity U upon action of an external force F 
1131 . 8 is dimensionless and compares the self-propulsion with drag- 
ging (some authors use the reciprocal efficiency 1 /S). The higher 8 
the more efficient the swimmer is. For an axisymmetric swimmer 
propelled along the symmetry axis, F- U = 1ZU 2 , where the scalar 
1Z is the appropriate hydrodynamic resistance. The work done by an 
arbitrary shaped swimmer and dissipated by viscosity in the fluid is 
given by 



V 



vdS = 2fi / E'.EdV 




Fig. 2. The swimming efficiency of the spheroidal microswimmer, <5, vs. the 
scaled elongation c/a (linear-log plot). The solid line corresponds to the exact 
calculation, the dashed line is the asymptotic result. 



where E is the rate-of-strain tensor, V is the fluid volume surround- 
ing the swimmer and S its surface. Expressing the product E:E 
as J]] uiiUJi + 2(diVj)(djVi) allows re-writing V for microswimmers 
self-propelled by purely tangential motions u as l 1 31 



V = n / wW + 2/i 



AS. 



[14] 



where k s = —(ds/ds)- n is the curvature measured along the path 
of the surface flow, expressible in terms of the unit tangential and 
normal vectors, s and n, respectively. Let us now estimate 8 of the 
spheroidal treadmiler described in the previous subsection. For a pro- 
late spheroid s = , n — e T , respectively, and k s can be calculated 
as 

1 dH ( 



H T Hr dr 



c(r 2_ C 2)3/2 



Since the solution 111 111 corresponds to irrotational flowQ, i.e. u> — 
the volume integral in 111411 drops out. Substituting the expression for 
the surface velocity II 1011 and k s into the surface integral in 111411 we 
find 



V = 4nfJ,CU 2 s(Ta - 1) {(1 + To) COth 1 T a - T a } 



[15] 



Collecting the expressions for the drag force Q5J, velocity of self- 
propulsion 111 111 and the dissipation 111511 one can compute the swim- 
ming efficiency, 



TZU 2 



_ 2{r a 2 -r a (r a 2 -l)coth- 1 r a } 2 

V ~ (r 2 -1){(1 + r 2 ) coth- 1 r a - T a f 



[16] 



8 is plotted as a function of the elongation c/a in Figure|2] Evidently, 
8 grows unbounded as c/a — > 00 does, corresponding, in the limit 
to a frictionless swimmer. For the spherical treadmiler 8 can be cal- 
culated from II14II with u = u s sin 9 eg, u s = |(7 and k s = l/o, 
V = 2 M f s § U 2 sin 2 9 i AS = 12 n^aU 2 . Dividing 6nfj.aU 2 
by V we find 8 — | in agreement with 111611 (see Figure |2]l and the 
theoretical bound (i.e. 8 < |) in 1131 . 

For the slender swimmer the asymptotic behavior of 8 can be esti- 
mated from II16II by expanding 8 in a series around T a = 1 and using 
where e 



l/VT 



1 



f «i. 



1 



2 (slogs) 



3 Since E 2 (ff m G T n J — the votricity ui is determined by the C n 



[17] 



> 4 terms in ^ 
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This result is shown in Figure[2]as a dashed line. For comparison, the 
efficiency of spherical squirmers self-propelled by propagating sur- 
face waves along their surface (the mathematical model of cianobac- 
teria fEl ) has the upper bound 5 < |, while numerically calcu- 
lated values of 5 do much worse than dragging and the corresponding 
swimming efficiency is usually less than 2% 1 13 1. Swimming by sur- 
face treadmilling is remarkably more efficient than the rotating helical 
flagellum 116], beating flexible filament 1 17], the Percell's "three-link 
swimmer" 1181 or locomotion by virtue of shape strokes 1 1411 1 91 . The 
surface treadmilling is probably superior to any inertialess swimming 
techniques proposed so far. 

Also, the swimming efficiency of the ellipsoidal treadmiler is 
superior by a factor of (clog 1/e) -1 over the estimate of 5 corre- 
sponding to the rod-like treadmiler with rounded ends derived from 
purely scaling arguments in the introduction. Therefore, the geome- 
try (via n s ) plays an important role in minimizing the dissipation in 
surface treadmilling, which is rather surprising since the drag force 
on slender nonmotile object does not depend on its shape to the first 
approximation. 

Optimal swimming. We can set an upper bound on 5 for a 
spheroidal microswimmer in terms of surface integrals of an arbi- 
trary velocity «(£) analogously to 1131 . The power dissipated in 
self-propulsion is bounded from below according to I1 1411 by 



-t--i 

V > 2{1 J U 2 K s dS = 4TVfJ,C T a (t, - l) J ■ 



2 (C) 



dc 



where we used the previously derived result for k s . The power ex- 
panded in dragging at the same speed is found from $9\\ and 111311 

as 



TZU 2 



2TVCIJ,T a 



(1 + r|) COth T a - T a 



^(OdC 



Combining the last two results we obtain an upper bound on 8 as 



5 < 



' u dC 



(7-2 - 1) ((1 + Tl) COth- 1 Ta -Ta) 



+ 1 

2 / : 



- dC 



The term in the figure brackets can be shown to be bounded 
from above by 2/3 while its maximum is obtained for u(() — 
Us yl — C 2 \f T a C 2 corresponding to the 2-term boundary velocity 
expansion [J^Jj with 62, &4 7^ , b m — 0, m > 6 (and, thus, repre- 
senting a rotational flow). Thus, for an elongated swimmer (r a — * 1) 
we arrive at 

<5 < 



1 



3e 2 log 1/e 

which does better than fl 1711 by a factor of (log 1/e) -1 and also su- 
perior over the scaling estimate for the rod-like swimmer by a factor 
of C(l/e). Note that the asymptotic behavior S ~ si ^ ^ji was 
derived from simple scaling arguments in the introduction. 

It can be demonstrated that the surface velocity II10II is not opti- 
mal, i.e. it does not minimize V for a prescribed propulsion speed U . 
To see this consider the slightly perturbed boundary velocity 



"(C) 



-2tw 



(t 2 -C 2 )MW 2 )- 



{U2G 2 (C)+M4G 4 (C)}, [18] 



such that I U4 1 <C \u2\ and U2 ~ u s . The solution of the linear prob- 
lem [|7l] yields |w| = 0(114) and therefore, the volume integral in 
ltT4ll is fj, J y u 2 AV = 0{u\). The surface integral in iTTHl can be 




( 

Fig. 3. The optimal boundary velocity u vs. the spheroidal coordinate C for 
elongation c = 2.5 a at various truncation levels : L = 4 (red) and L = 10 
(blue). The black line correspond to the one-term boundary velocity ftlOll 




c/a 

Fig. 4. The dissipation integral vs. the scaled elongation corresponding to 
optimal swimming at various truncation levels (linear-log plot): L = 4 (red) and 
L = 10 (blue). The black line refers to iflSTl - 



calculated as 

2/i / W 2 «: s dS' = U^I^Ta) + U2UaIa(To), 

Js 

where T^ija) is given by fl 1511 and T4 is some other function of 
T a . The velocity of propulsion can be found in the close form as 
U — U2J-2{to) + UiJ-4,(T a ) from [11311 . where T2 is equal to the ex- 
pression in the figure brackets in 111111 and JT 4 is some other function 
of T a . As the propulsion velocity to be fixed, we require U — u s J-2- 
This yields dissipation 

V = u s l 2 + u iUs (li - 2^ l 2 ^j + 0{u\). 

where the function in the brackets can be shown to be positive and 
bounded for all r a > 1, and vanishes only at r a — * 1. Therefore, 
one can always choose some u± < such that P < l2U 2 s leading to 
reduction in the dissipation in 111511 . The above perturbation analy- 
sis shows that, quite surprisingly, vorticity production could bring a 
reduction in viscous dissipation, leading to more efficient swimming. 

To address the question of optimal swimming we consider an ar- 
bitrary boundary velocity via the expansion, that meets all the above 
requirements for regularity and evenness in f , 



"(0 



-It a 



(r 2 -C 2 )Mi-C 2 )- 



u m G m (0- [19] 



where it follows from [JS1] that u m = —b m /(2c 2 T a ). To find a set 
of Fourier coefficients u m , m = 2, 4, 6, . . . L corresponding to the 
optimal swimming, one should minimize the dissipation integral, V, 
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while keeping the propulsion speed U fixed. The dissipation inte- 
gral V — — J s <r T ( u dS being bilinear in Uj, can be expressed as 
V = j^iii Ttj Uj . Note however that the tangential stress a T ( at 
the surface of the microswimmer requires the knowledge of the veloc- 
ity gradient at the surface (rather than velocity along). Alternatively, 
since the optimal velocity field is rotational, calculation of V from 
II 1411 requires the knowledge of vorticity everywhere. 

The propulsion velocity given by 111311 is linear in m, i.e. U = 
TjUj. The optimal set of coefficients m is to be determined from 

^-(■p — XU) = 0, or just from . Vij Uj = X Ti, where A is a 
Lagrange multiplier. We found the closed form optimal solution for 
the two- term boundary velocity [11911 . while for L > 4 closed form 
expressions are cumbersome, and numerical solutions were derived 
instead. Analogously to the theory for a 2-D swimmer (see the next 
section), where the explicit optimal solution was shown to acquire an 
infinite number of harmonics in the expansion for the boundary veloc- 
ity, increasing the truncation level L in II19II will further improve the 
efficiency of swimming, though the enhancement appears to be minor. 
To illustrate this, we calculate the optimal solution upon varying L. 
The optimal boundary velocity upon varying L is depicted in Figure 
[3] for the elongation of c = 2.5a and compared with the one-term 
expression II10II . The scaled dissipation integral, V/ fJ.au 2 , is plotted 
vs. c/a upon varying L in Figure [4] It can be readily seen that the 
convergence with respect to L is rather fast; the deviation between the 
results corresponding to L — 8 and 10 is less than 1% for all c/a and 
it vanishes at both limits c = and c/a — > oo. Thus, the 'intuitive' 



one-term boundary velocity f 1 1 Oil - that yields 5 ~ (eloge) - (see 
Figure[2]( is nearly optimal for a wide range of elongations and likely 
so for all elongations. 

2-D microswimmer. The two dimensional Stokes equations is con- 
veniently handled by employing complex variables 1 141 1 191 [2Tl 1221 . 
This allows explicit solution of the optimization problem for the el- 
liptical treadmiler. 

Denoting v — v x + iv y and d = \{d x — id y ), the Stokes equa- 
tions become 2fiddv = dp, Kedv = 0. The most general solution 
to this (with p real) is v = g-\- f — zg', p = — 4/iRe((?') where g, f 
are any pair of holomorphic functions 1201 . Solutions corresponding 
to multivalued g, f are also legitimate (provided the resulting v, p 
are single valued). It can be shown (using ||20ll below) that the mon- 
odromy of g (and of /) around a closed curve gives the total force 
exerted by the fluid on the interior of the curve. In particular in swim- 
ming problems this force must vanish and g, f are therefore always 
single valued. 

The element of force dF = dF x + idF y acting on a length ele- 
ment dz = dx + idy of the fluid can be expressed in terms of v, P 
and hence in terms of g, f. Straightforward calculation shows that the 
relation is 



dF = ip dz + (2iudv) dz = 2iu d(v - 2g). 



[20] 



Note that here (da;, dy) is tangent rather then the normal to the seg- 
menQ. 

We consider a 2-D swimmer shaped as an ellipse of semi-axes 
b,a — 1 ± a (with < a < 1) situated in the complex z = x + iy 
plane. It is then convenient to define a new complex coordinate £ by 
the relation z = £ + q/C- As £ ranges over the region \C\ > 1 the 
corresponding z ranges over the area outside the swimmer. In par- 
ticular the swimmer boundary corresponds to the unit circle £ = e l8 . 
Note that if we consider g, f as functions of £ rather then z then the 
general solution of the Stokes equations becomes 



1 - a/( 2 



[21] 



In the swimmer frame of reference, the boundary condition at 
infinity w(oo) = —U (where U is the laboratory-frame swimming 
speed) implies Laurent expansions 



f = J2 a ^~ n >9 = -u+J2 b ^~ 



[22] 



where U is arbitrarily appended to g. 

The boundary condition on the swimmer surface is fulfilled 
by matching to a prescribed boundary motion wL_ e is = 

w i®) = SnL-oa w ne m9 . It is useful to express w{6) as w = 

w_(C),C = e l8 where tu_ = 2^=0 w -nC™, w + = 
are both analytic outside the unit circle. Substituting 



w+(0 

1 122 II into 112 111 and matching on the unit circle we find 

C(i + «C 2 ) 



g(0 = w-(0, f(C)=w+(0 + 



c 2 



-(C)- 



In particular the swimming velocity is determined by the constant 
term in this expansion U = —too- The corresponding dissipation is 
calculated using 112011 as 

V = -RefMF = 2^fm(2 g - w) = ^ ± |n| M *. 



n— — oo 



Let us next focus on the case of an ellipse swimming by sur- 
face treadmilling. The boundary velocity w(9) being tangent to the 



swimmer boundary is expressible as w 



M0)=i(C-a/C)«(0) 



for some real-valued function u(9). Since we consider only swim- 
mers symmetric with respect to the x— axis, we assume u(6) to 
be an odd function allowing to write it as u = u n sin(n#) = 
i u n(C" — C~™)- I n terms of this the swimming velocity turn 
into U — —wo = |(1 + a)ui while the dissipation takes the form 

V = 2nn^2 n ((1 + a2 ) u n — 2«ii n _iu n+ i) . 

Which may also be written as V — | y~] VijUiUj for a corresponding 
tri diagonal matrix Vij. 

The optimal swimming technique for a given a is the one that 
minimizes the dissipation while keeping the swimming velocity U = 
I (l+a)wi fixed. The minimizer is the solution of -^-(P — Xui) — 
for Vi with A being a Lagrange multiplier, or just 

^^VijUj = X8i t i. [23] 

i 

It is readily seen that the coefficients with even k are not rele- 
vant to the optimal swimming and should be set to zero to minimize 
viscous dissipation. (This is also clear from the fact that U2k corre- 
spond to flows which are antisymmetric with respect to the y— axis.) 

Denoting = lt2fc+i, fc > and writing A = 2u s 7r/i(l + a 2 ) 
with u s an arbitrary normalization constant having dimensions of ve- 
locity we obtain from 112311 the recursion relation 

(2k + l)b k - £(&& fc _i + (k + l)b k+1 ) = u 3 8 kt0 , 



where £ 



Multiplying by x k and summing over k this 



transforms into a differential equation for the generating function 

1 , 



+ ^h'{x)B{x) + 1 



0. 



h{x)B'{x 

where we defined h(x) = 
to this is B(x) = 



/h(x) JL - y/h(x') 

tegration. Requiring the coefficients b k to decay for large k implies 




a) (a; — —). The general solution 
where C is a constant of in- 



4 The sign convention here is that dF is the force exerted by the l.h.s of the (oriented) segment 
dz on its r.h.s. 
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that B(x) must be analytic inside the unit disc and hence its potential 
singularity at x — a must be avoided. This determines the integration 
constant to be C = a, so that we may write 



B(x) = 



1 



Ax' 



log 



/xZ + 



where x+ — — — x and X- — a — x. The corresponding swimming 
velocity and the dissipation are, respectively, 



u = <i± ^il = £(i + a)B ( ) = log ( i±£ 

2 2 At; \ 1 — a 

V = — Pij U i U 3 = \^ Ul = n f M (^ + « 2 )"sS(0) = 



— u s (l + a log 
2a \ 1 — a 

Therefore, combining the last two results yields 

2a 47r/i[/ 2 



V 



(l + «) 2 lo f 



(^) 



We recall that in 2-D the dragging problem admit no regular solution 
within the Stokes approximation Thus defining the swimming effi- 
ciency as S = (F ■ U)/P = 1ZU 2 /V makes no sense in the present 
2-D context in which F, 1Z are not defined. This may be considered 
as a mere issue of normalization. We therefore use here an alternative 
definition of swimming efficiency 1 19 1 where 



47rp[/ 2 _ (1 + a) 



V 



2a 



■ log 



1 - a 



[24] 



In the slender limit, a — > 1, or, % = — e 



0, as the el- 



_ i- 

b ~ 1 + a 

lipse degenerates into a needle, the efficiency grows logarithmically 
unbounded as 

<r~2io g (i/ s ). 

It may be of interest to note that truncating our expansion to include 
any finite number of modes would lead to B(x) which is not only 
polynomial in x but also algebraic in a. This then implies that the 
(truncated) efficiency S* oc B(0) would be algebraic in a implying 



that 112411 must be modified to a bounded expression. Thus (in contrast 
to the 3-D case) one cannot obtain the correct asymptotic efficiency 
without retaining all the modes. 

The optimal boundary velocity may be found explicitly as 



w{6) 



dz 



]T u k sin(fc0) = ^| bk 



(2fc+i)e 



i(2fc+l)S\ 



dO 

dz 



2i 



u s ^Im{e ie B (e 2ie )}. 



Using the explicit expression we have for B(x) and 4| we find the 
absolute value of w is given by 



1 + a 2 , Va - c 2i » + JTTc 

Us j=— log 



v/l/a — a 

while its direction is known to be tangential. Thus, in the limit e — *• 
one finds (provided \tv — 8\) that 



\w/U\ 



1 + 



log 2 sin G\ 
logl/£ 



and the boundary velocity approaches the constant U though only at 
a logarithmic rate. 

Concluding remarks. In this paper we examined the propulsion of 
elongated microswimmer by virtue of the continuous surface tread- 
milling. As the slenderness increases, the hydrodynamic disturbance 
created by the surface motion diminishes, i.e. the microbot is pro- 
pelled forward with the velocity of the surface treadmilling, while 
surface, except the near vicinity of the poles, remains stationary in the 
laboratory frame. As a result of that, the 'cigar-shaped' treadmiler 
is self-propelled throughout almost quiescent fluid yielding very low 
viscous dissipation. The calculation of optimal hydrodynamic effi- 
ciency of the 3-D and the 2-D microswimmers reveals that the pro- 
posed swimming technique is not only superior to various motility 
mechanisms considered in the past, but also perform much better than 
dragging under the action of an external force. 
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